Spectra of massive and massless QCD Dirac operators: A novel link 
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We show that integrable structure of chiral random matrix models incorporating global symme- 
tries of QCD Dirac operators (labeled by the Dyson index /3 = 1,2, and 4) leads to emergence of 
a connection relation between the spectral statistics of massive and massless Dirac operators. This 
novel link established for /3-fold degenerate massive fermions is used to explicitly derive (and prove 
the random matrix universality of) statistics of low-lying eigenvalues of QCD Dirac operators in 
the presence of SU(2) massive fermions in the fundamental representation {fi = 1) and SXJ(N C > 2) 
massive adjoint fermions (f3 = 4). Comparison with available lattice data for SU(2) dynamical 
staggered fermions reveals a good agreement. 
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Explicit knowledge of spectral statistics of low-lying 
eigenvalues of the Dirac operator is required to un- 
derstand the phenomenon of chiral symmetry breaking 
(xSB) in quantum chromodynamics (QCD). It has first 
been conjectured by Verbaarschot and collaborators 
that extreme infrared limit of the QCD Dirac opera- 
tor spectrum can be described by the large-iV chiral 
Random Matrix Theory (RMT) that models the true 
Dirac operator T> by N x N block offdiagonal matrix 
£>RMT _ ffdiag(iW, iW^), W being n x m rectangular 
random matrix [see Eq. (Q) below]. In such a formula- 
tion, N = n+m is an analog of dimensionless space-time 
volume V, while v = \n — m\ is equivalent to the topo- 
logical charge (equal to the number of zero modes of T>) . 
If, in addition, the entries of W are chosen to be real, 
complex, or quaternion real, the random matrix £> RMT 
possesses proper antiunitary symmetry (labeled by the 
Dyson index f3 — 1,2, or 4) and, hence, correctly re- 
produces both the underlying symmetries of the Dirac 
operator and the x^B pattern associated with it. On the 
language of the chiral QCD Lagrangian, the above ap- 
proach corresponds to the limit [| 1/A< y 1 / 4 <c l/m^ 
(A is a typical hadronic scale and is the pion mass) in 
which the kinetic term in Lagrangian can be neglected, 
and only the global symmetries of the Dirac operator 
become important. Recently, RMT phenomenology has 
been put onto a firm field theoretic ground represented 
by the framework of finite-volume partition functions jjj 
and the partially quenched perturbation theory Q. 

On the microscopic scale ~ 1/VE, chiral RMT (defined 
for a given topological charge v M) leads to parameter- 
free predictions for the unfolded microscopic spectral 
density p s (X) = limv^^iVY,)- 1 p(\/VT,) of the Dirac 
operator. Here, the absolute value E of the chiral con- 
densate (the order parameter of x$B) is related to the 
Dirac spectral density p(0) at zero virtuality through the 
Banks-Casher relation E = Trp(0)/V Q. In a series of 
papers , RMT predictions have been confronted to the 
lattice data, and good agreement has been found for the 
spectral density, two-level correlation function, and dis- 



tribution of the smallest Dirac eigenvalue for massless 
lattice data of all xSB patterns. The lattice data for 
massive fermions have recently appeared B as well. 

Unfortunately, available theoretical results for micro- 
scopic spectral correlators of massive QCD Dirac oper- 
ators are rather poor, being restricted to the (3 = 2 
symmetry class |9|,|l0| associated with the gauge group 
SU(A C > 3) in the fundamental representation. It is 
the aim of the present Letter to show that integrable 
structure of chiral RMT results in a simple but powerful 
link between the spectral statistics of massive and mass- 
less QCD Dirac operators for all three symmetry classes 
f3 = 1,2, and 4. The connection relation [Eq. (g)], es- 
tablished below for /3-fold degenerate massive fermions 
within the framework of chiral RMT, relates partially 
unknown massive spectral correlation functions to the 
massless ones (taken at both positive and fictitious nega- 
tive energies) pT| • As the latters have already received a 
detailed study in the literature, this link not only solves 
the problem posed but also provides a particularly simple 
proof of RMT-universality of massive correlation func- 
tions, that becomes a consequence of celebrated univer- 
sality [iJfQ 

proven for the massless case. 
Let us start with the definitions (l^Jl^]. The joint 
probability distribution function of chiral random matrix 
ensemble associated with Nf massive fermions in the sec- 
tor of topological charge v is given by 
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{A} > are the eigenvalues of the matrix 
A„ ({A}) = n"<j(^j — ^j) 1S ^e Vandermonde 
determinant, the weight function wp^ is wp.v(X) = 



Here. 
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, V(X) is the finite-polynomial confine- 
ment potential, and the topological charge v is taken to 
be positive integer or zero. 
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The p-point correlation function in the above ensemble 
is expressed as [|16| 



R 



(Ai,...,A p ) — 



(n — p)\ 
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, X p , X p+ i, . . . , X n ). 



For p — this yields the mass-dependent partition func- 
appearing in Eq. ([!]). The unfolded 
spectra of the Dirac operator are then obtained from the 



appropriately unfolded spectra R 



({A}) of associ- 



ated random matrix model Eq. @j) by a simple transfor- 
mation of variables: 



ps(X 1 ,...,X p ) = 2^H\X l 



\Ri Nf ' u ' m (Xl. 



A 2 ). (3) 



k=l 



We notice that for massive correlation functions this also 
demands to rescale the quark masses, /x/ = m/V£. This 
completes our definition of the model. 

In what follows, we assume that the massive fermions 
are /3-fold degenerate, M.p — {m^p, . . . , m,N f 1/3 ), so 
that appropriate matrix ensemble is given by the joint 
probability distribution function Pn Nf ' ]/ '^ ({A}). Since 



({A}) contains a positive definite 
)P , it can conveniently be absorbed 



in this case, Pn 

factor Il/^i^ + TO / 
into a Vandermonde determinant of a larger dimension 

An+Ar, ({A}, {-to 2 }) = A„ ({A}) A Nf ({-to 2 }) 

n N f 

xIHI( A ' + m /)- ( 4 ) 
i=i/=i 

This immediately results in a pretty fact that the par- 
tition function of the model Eq. (Q) with (3Nf massive 
fermions can be expressed in terms of associated massless 
Nf- point correlation function R^+n) N (i~ m2 }) °f the 
matrix ensemble of larger dimension, (n + Nf) x (n + Nf), 
taken at fictitious negative energies: 
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The same strategy is applied to the p-point correlator, 
Eq. (0). After a few transformations, we arrive at the 
following remarkable relationship: 
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Finite-™ Eq. (|S|) establishes a link between massive and 
massless spectral correlators via associated chiral ran- 
dom matrix ensemble, and represents a basic relation to 
be examined in the rest of the Letter, where we consider 
the most interesting symmetry classes /3 = 1 and 4. It 
should be stressed that, in spite of seeming simplicity, 
the link Eq. (fj) is not obvious as it involves correlation 
functions of massless chiral ensemble taken at both pos- 
itive and fictitious negative energies. We also wish to 
emphasize that, in the microscopic limit, Eq. (0) imme- 
diately leads to the RMT-universality of the microscopic 
massive correlators which becomes a simple consequence 
of the universality phenomenon established for massless 
correlation functions [[l2]jl]| . 

(i) Let us turn to the (3 = 4 symmetry class, associ- 
ated with SU(A C > 2) massive adjoint fermions. In the 
vicinity of the hard edge, the unfolded p-point spectral 
correlators in the massless chiral model Eq. (p i), N f = 0, 
admit quaternion determinant representation ]]! q| 

4°' M) ({A}) - Qdet [/ 4 (A l; A,)]^ ^ . (7) 

Here, the 2x2 matrix kernel f± = fp—4 [§~7j 



M ' j ~ v. lp(X,Y) S P (Y,X) J 



Dp(X,Y) = -d Y S (3 (X,Y), Ip(X,Y)= dZS (Z,Y)-e(X-Y)5p A , (8) 



x 



e(X) = (l/2)sgn(A), is uniquely specified by the func- 
tion |il 



S 4 (X,Y) = 2K 2v+1 (2X x l 2 ,2Y x ' 2 
J 2 „{2X^) f 2Yl/2 



, Yi , L dtJ 2u+2 {t) (9) 

with K a {X, Y) being the Bessel kernel @: 

is (v v\ x Ja+i{X)J a {Y) -YJ a+ i(Y)J a (X) 

K«(X,Y) = 2 (X2-y 2) • ( 10 ) 



It is important to stress that Eqs. (|^) and ([To]) hold 
generically [p"3| for arbitrary finite-polynomial confine- 
ment potential V(X) in Eq. (|l|). 

As long as at j3 = 4 the definitions Eqs. ([!]) and (||) 
do not discriminate between the positive and (fictitious) 
negative eigenvalues Afe, Eqs. @ and (|h]) remain valid 
at negative arguments also, provided one makes use of 
the fact that J a (iX) = i a I a (X). This circumstance al- 
lows us, at once, to derive closed expressions for p-point 
massive correlation functions in the microscopic limit. 
Straightforward calculations based on Eqs. @, (||) and 
(0) yield: 
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FIG. 1. Microscopic massive density p s (A; /il 



2 4 6 8 '0 12 A '"0 2 4 6 8 '0 12 A 
Eq. (12), plotted versus lattice data simulated in Ref. [8] for 8 flavors of 



masses n = 8.74, 6.31, and 4.26. For discussion of species doubling and systematic deviations at large A we refer to Ref. [8]. 
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/ 4 (Af,Aj) h(Xi-nj.) 



(11) 



Qdet[/ 4 (-^,-4)] 

1 < *i j < P, 1 < /■ /' < ^V/- We remind that this result 
applies to 4-fold degenerate quark masses, Uf = rrifVT,. 



Microscopic density for the particular case of 4 de- 
generate fermions of mass /1 is of special interest as it 
can be compared to available lattice data for dynamical 
SU(2) staggered fermions in the fundamental representa- 
tion simulated in Ref. Q at v = 0. [Because of the lattice 
symmetry of staggered fermions they belong to the sym- 
metry class (3 — 4:]. Computing quaternion determinants 
in Eq. (0), we come down to 



PsWl4) 



2|A| S 4 (A^,A 2 ) 



S 4 (-,j 2 , A 2 )5 4 (A 2 , —fi 2 ) - I 4 (A 2 , - A t 2 ) J D 4 (A 2 



(12) 



We have explicitly checked that the limit fi — > repro- 
duces the known massless result [Q at a shifted topolog- 
ical charge v — ► v + 2 [^0| . Theoretical results plotted 
in Fig. 1 for three different values of /1 show reasonable 
agreement with numerical data. 

We close our consideration of (3 = 4 symmetry class by 
giving a compact expression for the previously unknown 
massive RMT (or finite-volume pl"l) partition function 
with 4-fold degenerate massive fermions [see Eq. (B)]: 



, /xat / 1 4 ) 



(-l) jV /Qdet[/ 4 (-/i 2 ,-/x 2 / ,)] 



iA W/ (-^)i 4 n&^ 



(13) 



Here, only nontrivial mass dependence has been dis- 
played. 

(ii) Now we turn to the /3 = 1 symmetry class asso- 
ciated with SU(2) massive fermions in the fundamental 
representation. In this case, the modulus of the Vander- 
monde determinant in Eqs. (|j) and (||) makes allp-point 
correlation functions to be nonanalytic functions of their 
arguments. This is exactly the reason of why one cannot 
use known expressions for massless correlation functions 
^'"'^({X}) to naively compute them at negative ener- 
gies. Below we show how to circumvent this obstacle for 
the simplest situation of the spectral density with a single 
quark mass. Extension to higher order correlation func- 
tions and/or larger number of masses is straightforward. 

In accordance with the connection relation Eq. (|^), 
the finite-n massive spectral density equals 



R^f\x,- m 2 ) 



R 



(0,^,1) 



( — TO 2 ) 



(14) 



where, for definiteness, we have fixed the dimension of the 
massless random matrix ensemble to be even, 2n; from 
now on, the superscripts are omitted for brevity. We ob- 
serve that the function i?2n,2(A, —to 2 ) can be evaluated 
through the functional derivative of i?2n,i(— m 2 ', [W]) = 
R2n,i(—m- 2 ) with respect to the confinement potential 
W, exp{-W(X)} = A(' y - 1 )/ 2 exp{-l/(A)}: 



Ri 



2 (A,-m 2 ) = i? 2 „,i(-m 2 ) (i? 2 „,i(A) -5{X + m 2 )) 

- ww R ^ { - m2 ' [w]) - (15) 

To facilitate taking the functional derivative in Eq. (|lj) , 
we utilize the approach of Ref. |^2| [see Eq. (A6) of sec- 
ond reference], but express i?2n,i(— ™ 2 ; [W]) in terms of 
arbitrary polynomials Pjix) rather than in terms of the 
skew orthogonal ones: 



R2nA~rn 2 ;[W}) = -e 



-w(- 



2n-l 

2 ) ]T P] (-m 2 )n ]k [W] 

j,k=0 



x / dZ e- w ^ Pk (Z). 



(16) 



The In x 2n real antisymmetric matrix fij k [W] is the 
inverse to the matrix 17 



M 



-\-oc 



dxdy e~ w ^' w ^e(x - y) Pj (x)p k (y), (17) 
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e(x) = (l/2)sgn(x). Substituting Eq. @ into Eq. tffify, 
and then into Eq. we are able to express the finite-n 

massive spectral density i^n— 1 lW m tne f° rm 



s[ 2n) (\,\) 



S[ 2n \-m 2 ,\)s[ 2n) (\,0) 



(18) 

7f n) (0,A)£)f n) (-m 2 ,A) 



sf n) (-m 2 ,0) 
that contains the entries of finite-n, 2 

(2n), 



ff=\{X,Y) of the mass- 
less ensemble [see Eq. (g)]. In deriving Eq. 
(O) we have used both the representation |T^] 

Z)e- w ( z 1pk(Z), and Eq. (§). Finally, taking into ac- 
count Eqs. (||), and the universal jl3) formula Q 



trix kernel ff n '(X,Y) 



S 1 {X,Y) = K v _ l {X 1 / 2 ,Y 1 / 2 ) 

MX 1 / 2 ) ( f Yl/Q 



dtJ v _ 2 {t) 



(19) 



(valid in the vicinity of the hard edge) , we deduce a closed 
expression for the microscopic single-mass spectral den- 
sity pg (Aj/u). It exhibits the quaternion determinant 
structure of Eq. ([l2|) with obvious changes in arguments 
of 5, D, and I functions as is given by Eq. (fl8]). As 
a consistency check, we have verified that for fi = it 
reduces to the known result [EJ for one massless flavor. 
Let us stress, that universal form |b| of the function 
Si thus confirms the universality of the massive spectral 
density following in a more general context directly from 
the microscopic limit of the connection relation Eq. (||). 

In conclusion, we have derived universal expressions for 
spectral correlators of massive chiral matrix ensembles 
corresponding to /3-fold degenerate massive fermions, by 
establishing a new link between the statistics of mas- 
sive and massless random matrices. The results obtained 
have been compared to the available lattice data associ- 
ated with (3 = 4 xSB pattern in low-energy QCD. 

We thank the authors of Ref. || for making their data 
available to us. E.K. acknowledges early discussion with 
S.M. Nishigaki. The work of G.A. was partially sup- 
ported by EU-TMR grant No. ERBFMRXCT97-0122. 

Note added.-After completing this work, the preprint 
p3| by T. Nagao and S.M. Nishigaki on finite-volume 
partition functions has appeared. In particular, these 
authors give alternative representations of our Eqs. jll| ) 
and (|J). 
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